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Abstract 

An extension of the adiabatic factorization of the time evolution operator is studied for spin in a 
general time varying magnetic field B(t). When B(t) changes adiabatically, such a factorization reduces 
to the product of the geometric operator which embodies the Berry phase phenomenon and a usual 
dynamical operator. For a general time variation of B(t), there should be another operator N(t) in 
the factorization that is related to non-adiabatic transitions. A simple and explicit expression for the 
instantaneous angular velocity of this operator is derived. This is done in a way that is independent of 
any specific representation of spin. Two classes of simple conditions are given under which the operator 
N(t) can be made explicit. As a special case, a generalization of the traditional magnetic resonance 
condition is pointed out. 



1 Introduction 



The purpose of this paper is to study a unique factorization of the time evolution operator for spin in a 
general time varying magnetic field that should be seen as a generalization of the adiabatic factorization. 

Let us first recall that in the proof of the quantum adiabatic theorem as presented in standard texts such 
as Messiah [I], a path-dependent operator A(t) and a dynamical operator D(t) are successively extracted 
from the time evolution operator. Then, there is left an operator N(t), the product of which with A(t) and 
D{t) should equal the time evolution operator U(t), i.e., U(t) = A(t)D(t)N(t). The rest of the proof of 
the quantum adiabatic theorem amounts to the proof of the statement that N(t) approaches the identity 
operator in the adiabatic limit. The nature of the path-dependent operator has been clarified with the 
concept of the geometric phase [2J El H] - Near the adiabatic limit, the operator N(t) (we shall sometimes 
refer to it as the non-adiabatic operator) can be analyzed in the spirit of perturbative analysis 5, 6, 7 . For 
a general time variation of the parameters in the Hamiltonian, the non-adiabatic operator N(t) in the time 
evolution operator is responsible for transitions among instantaneous eigenstates. While the operator A(t) 
is characterized by the geometry of the path traversed by the set of parameters in the Hamiltonian, and 
D{t) is the usual dynamical operator, a simple characterization of the operator N(t) that corresponds to a 
general time variation of the parameters in the Hamiltonian seems to be still lacking. 

This paper intends to clarify the structure of the operator N(t) for spin in a general time varying magnetic 
field from a perspective that is independent of perturbation theory. The question we ask is how the geometry 
of the path of the direction of the magnetic field and the dynamical effect contribute in determining N(t). 
As it turns out, N{t) in this case is determined by an angle parameter associated with the geometry of 
the path, a dynamical angle that characterizes D(t) and the speed function that describes how the path is 
traversed in time. While this is reminiscent of the fact that the operator N(t), which is an SU(2) rotation, 
should in principle be characterized by three parameters, our point in this paper is that the instantaneous 
angular velocity of N (t) can be explicitly expressed in a simple way by using the time dependent parameters 
associated with the Hamiltonian. The factorization of the time evolution operator U (t) derived here can 
be applied to the most general time variation of the magnetic field, including the adiabatic scenario as a 
special case. Through this natural factorization of U(t), we find that there exist new explicit solutions of 
time evolution that go beyond the scope of perturbative analysis. 

In a certain sense, the approach we adopt in this paper differs from the traditional adiabatic perturbation 
theory (See, for example, [6] for a review). Let us observe that in the proof of the adiabatic theorem pQ, the 
path-dependent operator is formally constructed out of the energy eigenstates, i.e., 



Such a formal expression in terms of eigenstates is necessary because in the general situation no specific 
information of the Hamiltonian is available. The traditional adiabatic perturbation theory is strongly in- 
fluenced by the use of instantaneous eigenststes. While such an approach is quite general in nature and 
have wide applications, it can also lead to complicated calculations and the results obtained from such an 
approach cannot be extrapolated to situations that are not within the scope of perturbation theory. We are 
of the point of view that in analyzing a specific problem such as the spin in a time varying magnetic field, 
the use of energy eigenstates can be avoided and that the time evolution can be studied by exploring the 
algebraic structure of the Hamiltonian alone in constructing a factorization of the time evolution that can 
then be applied to any representation. 

The Hamiltonian for a spin S with magnetic moment — fcS in a time varying magnetic field is 




(1) 



where \ip m (t)) is the instantaneous eigenstate of H(t) that satisfies 



(tp m {t)\iJm(t)) =0. 



(2) 



H = kB(t) ■ S = kB{t)n{t) ■ S, 



(3) 



2 



where A; is a constant, B(t) = B(t)n(t) is the magnetic field and n(i) is a unit vector. 
If B(<) changes adiabatically, the time evolution operator takes the form 

U A {t) = A(t)D(t), (4) 

where ^ 

A(t) = Texp(-i I (nxnj-Sdr), (5) 
Jo 

and 

D(t) = exp ( - i I fcS(r)n(0) • S dr) . (6) 
Jo 

The operator A(t) is determined by the path traversed by the direction of the magnetic field n(i) on the 
unit two-dimensional sphere. (We set h = 1.) 

For the spin 1/2 case, the operator A(t) is written in this way using Pauli matrices in Berry's iterative 
approach to calculate non-adiabatic corrections to the geometric phase [5]. In the factorization of time 
evolution studied in the present paper, the two factors A(t) and D(t), and another factor N(t) to be derived 
later belong to an arbitrary finite-dimensional irreducible representation of S£/(2), so the method and result 
here apply to any spin. 

First, let us fix certain notations. Recall from the quantum mechanics of spin that the dot product in the 
Hamiltonian indicates an isomorphism between the three dimensional vector space and the Lie algebra with 
a certain basis: {— iSi, —iS2, —iSs}. To be more specific, { — iSi, —1S2, —iSs) here is a basis of an arbitrary 
finite dimensional irreducible representation of the Lie algebra su(2). So, in effect, the operators that are 
concerned in this paper can be thought of as finite dimensional square matrices. 

Let ei(t), e2(t), and es(t) be three orthogonal unit vectors determined by the differential equations 

Bi{t) = (n x n) x Bi{t), (7) 

with i = 1,2,3, and the initial condition 

ei (0) = n'(0), e 2 (0) = n(0) xn'(O), e 3 (0) = n(0), (8) 

where the prime in n'(0) means derivative taken with respect to arc length. The vectors ei(t) and e2(t) 
satisfy the condition ei(i) ■ e2(t) = 0. When viewed as tangent vectors, they have the geometric meaning of 
being parallel transported along the curve n(i) on the unit two-sphere. 
Let us define Si, S 2 & n d S3 as follows: 

Si=ei(0)-S, S 2 = e 2 (0)-S, S 3 = e 3 (0) ■ S. (9) 

In accordance with the quantum mechanics for spin, these operators are assumed to satisfy the usual com- 
mutation relations [Sfc,S;] = ieklmSm- 

Finally, the operator A(t) is, by definition, the solution to the differential equation 

A(t) = H(n x h) • S)A(t), (10) 

with the initial condition that ^4(0) is equal to the identity. Because A~ X A + A~ x A = 0, we have A = — 
A~ 1 AA~ 1 . Then it follows from the definition of A(t) that 

A- 1 {t) = A- 1 (t)(i{nxh)-S). (11) 

2 Mathematical Framework 

In this section, we give an introduction to existing knowledge on the fundamental properties of the operator 
A(t) and also a rotation angle between two frames naturally associated with the path traversed by the 
magnetic field. This preparation clarifies the general factorization of time evolution completed in the next 
section. 
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2.1 Two properties of A(t) 

It is well known [3] that the Berry phase has mathematical interpretations in terms of fundamental concepts 
in differential geometry and Lie groups. Here we review the relevant mathematics from the point of view 
of the operator A(t). As we shall see, A(t) plays an important role in deriving the general factorization of 
U{t). The two characterizing properties of A(t) are as follows. 
Property 1. First, we have the following property of A(t): 

A- l {t)ei{t)-SA{t) = e 4 (0)-S = 5 j; . (12) 

To prove this, note that this clearly holds for t ~ by definition of A(t). Now it suffices to prove that the 
left hand side is constant, or, its time derivative is zero. This can be verified by using the expressions for 
A{t) and A _1 (i), the equation ei(t) = (n x h) x e,(t), and the formula [a • S,b • S] = i(a x b) • S. 

This property of A{t) [5] is best understood in terms of the adjoint representation of the Lie group 
G generated (through the exponential map) by the Lie algebra with the basis {— iS\, — 1S2, — iSa}. So G is 
isomorphic to SU(2) or 50(3). In the case where G is isomorphic to SU(2), it is well known that this adjoint 
representation is a 2 to 1 covering map from G to SO(3), and so in general the holonomy associated with 
the standard metric on the two sphere is insufficient in determining the Berry phase information contained 
in A(t). So the following property is indispensable if one wishes to cast the solid angle Berry phase from the 
point of view of A(t). 

Property 2. For a cyclic change of n(i) with n(T) = n(0), we have A(T) = cxp[— iiln(0) ■ S] , where fi is 
the solid angle subtended by the closed path at the origin of the parameter space. Here the parameter is B. 

This is equivalent to associating a solid angle Berry phase to each individual energy eigenstate. To prove 
this property, let \ip m (n(0))) be an initial eigenstate of n(0) • S, i.e., 

n(0) • S|V>m(n(0))) = m|V> ro (n(0))>. (13) 

Then it can be shown, using Property 1, that \tp m (n(t))) = A(t)\ilj m (n(0))) satisfies 

n(t) • S|^ m (n(t))> = m|V m (n(t))), (^ m (n(t))|^ m (n(i))) = 0. (14) 

Now one may choose a locally single valued |m(n)) such that \ip m (n(t))) — e 1 " 1 ^ |m(n(t))). Then e* 7 ^ for a 
closed path can be calculated to be e ~ lm,a by Berry's phase two-form method. Since the |^> m (n(0)))'s form 
a complete basis for the representation space, it is clear that A(T) = exp[— ifin(O) • S]. 

In conclusion, although the operator A(t) has been written only as a time ordered exponential, it never- 
theless should be seen as an operator that can be explicitly constructed for a given representation of SU(2) 
and for a given n(t). This can be seen as a consequence of the fact that the instantaneous angular velocity 
has the special form (n x ri), so that the above two properties hold for A(t). In a specific representation, 
one has A(t) = J2 m IV'm(i))(V'm(0)|. Because the matrices S\,S2, and S3 can now be written explicitly, one 
can obtain explicit solutions to the eigenvalue equation n(i) • S|^ m (n(i))) = m|^ m (n(t))). The condition 
(tp m (n(t))\4'm(n(t))) = then fixes a phase factor for the eigenfuction and A(t) can then be made explicit. 
This contrasts the non-adiabatic operator N(t) to be discussed later: N(t) can be made explicit only in 
special cases, although the instantaneous angular velocity of N(t) has an explicit expression. 

2.2 Relative rotation between two triads 

In using an iterative method to calculate non-adiabatic corrections to the spin 1/2 geometric phase, Berry 
considered the relative rotation between the two triads (ei(t), e2(t), ^(t)) and (n'(t),n(t) x n'(t), e3(i)) [5]. 
The relation between the two triads can be expressed as 

ei(*) = n'(t)cos/3 + n(i) x n'(i) sin/3, (15) 

e 2 (f) = -n'(i)sin/3 + n(£) x n'(i)cosj8. (16) 
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The condition e'i ■ e 2 = then implies that 

P'(t) = -(n(t)xn'(t)).n"(t). (17) 

So /3(t) can be written as 

Mt) 

(3(t) = / -n" • (n x n')ds, (18) 
Jo 

where l(t) is the arc length traveled through by n(t). Thus, (3(t) is a geometric quantity in the sense that 
it is determined by the path traversed by n(t) on the two-sphere and is independent of how the path is 
parameterized. 

If l(t) increases with t, then the relation between n'(t) and the time derivative h(t) is h(t) = n'(t)\h(t)\. 
Then we also have 

/9(t) = / o ' ( n x ( 19 ) 

•>° | n | 

So /? is explicitly determined if n(i) is given. 

3 General factorization of time evolution 

In this section we study how the adiabatic factorization can be generalized when the magnetic field varies in 
a general way. The purpose is to derive an explicit expression of the angular velocity of the non-adiabatic 
operator in terms of the relative rotation f3(t), which is a geometric quantity and the dynamical quantity 
4>(t) = — J kB(r) dr in a way that is independent of specific representations of spin. 

Let the non-adiabatic operator be N(t) so that the time evolution operator in the general situation is 

U{t) = A{t)D{t)N{t). (20) 

Substituting this expression into the Schrodinger equation 

iU(t) = H(t)U(t) (21) 

we obtain 

iA~ 1 ADN + iDN + iDN = A^ 1 HADN. (22) 

By the first property of the operator A{t) from last section, it is clear that iDN = A~ 1 HADN, so that we 
have 

A- X ADN + DN = 0, (23) 

or, 

N= - D- 1 A- 1 ADN = iD- 1 A- 1 {n x n) • SADN. (24) 

This equation determines the operator N(t) that satisfies the initial condition N(0) = I. 
To calculate iD~ 1 A~ 1 (n x ft) • SAD in the above expression, let 

n(t) x n(t) = a(t)ei(t) + b(t)e 2 (t). (25) 

Then the relative orientation of the two triads from last section tells us that 

a(t) = |n(i)|sin/3(t), b(t) = \n(t)\ cos/3(i). (26) 

By Property 1 for A(t), we have 

A^A= -i(a(t)S 1 + b(t)S 2 ). (27) 

Using the formulas 

e -i4>S 3 g ie i4>S 3 = Si C0S(j) + g 2 sin(j) ^ (28) 
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e -#S 3 £ 2e #S 3 = 5 1 (_ sill 0) + S 2 coscj), (29) 

we now have 

- D~ l A~ x AD = ia{t){S x cos<j) + S 2 sin^) +ib(t)(Si(- sin (f>) + S 2 cos^>). (30) 

where ^ 

cj)(t) = - I kB(r) dr. (31) 
Jo 

The expression for — D~ x A~ x AD can now be written in the form 

- D- X A~ X AD = i5i|n(t)| sin[0(i) - #(t)] + iS 2 \n(t)\ cos^(t) - (32) 

from which we get 

7V(t) = Texp(-i^ (wi(r)Si + w 2 (r)5 2 )dr) , (33) 

where the two non- vanishing components of the instantaneous angular velocity are 

wi(t) = -|n(f)| sin[/3(<) - 0(f)], w 2 (f) = -|n(f)| cos[/?(f) - 4>{t)}. (34) 

The expression for N(t) in terms of the explicitly determined angular velocity is the main result of this paper. 
|ri(t)| here has the simple meaning of being the velocity at which the path n on the two-sphere is traversed 
in time. Together with |n(t)|, the two naturally defined quantities /3(f) and 4>(t), one geometrical and the 
other dynamical, provide a simple expression for the angular velocity which in turn completely determines 
the operator N(t). 

Since Si and S 2 can be expressed in terms of the raising and lowering operators S+ = Si + iS 2 , and 
S- = Si — iS 2 , N(t) also can be written as 



7V(<) = Texp(^ (z(T)S+-z*(T)S-)dT), 



(35) 



where 



z{t ) = Mi e *Ifl(*)-*(*)]. (36) 

In a certain sense, N(t) here is analogous to the non-adiabatic operator that we derived for the Landau 
problem with a general time dependent electric field, where it is a time-ordered coherent state displacement 
operator that describes inter-Landau level transitions [5]. However, the time-ordered coherent state displace- 
ment operator there can be explicitly determined because the commutation relation [a, a'] = 1 is relatively 
simple. In fact, N(t) for this time varying Landau problem can be shown to be equal to a numerical phase 
factor times the ordinary exponential. In contrast, the time-ordered exponential here has no immediate 
relation to the ordinary exponential. This is due to the fact that we have a different commutation relation 
between S+ and S-, i.e., = 2S3. We shall demonstrate below that N(t) can be made explicit for 

special cases where additional conditions are specified. 



4 Two classes of explicit factorizations 

The general factorization includes previous studied cases such as the adiabatic approximation and the sudden 
approximation as special cases. It may also be applied to study in detail certain generalizations of the 
quantum adiabatic theorem such as when the parameter passes through the degeneracy point [TU] ■ 

Let us take the sudden approximation for example. It states that if the parameters in the Hamiltonian 
changes (even if appreciably) in a sufficiently short time, then under quite general conditions, the quantum 
state remains close to the initial state pQ. Suppose now that the magnetic field keeps its magnitude B 



G 



constant but changes direction in a very short time, then the dynamical quantity (p(t) = — J* kB(r) dr is 
close to zero (because t — > 0) and thus can be ignored. We thus have D{t) = 1. If <f)(t) = 0, we have 

N(t) = Texp(-i^ (co 1 ( T )S 1 +uj2{T)S2)dT S j, (37) 

= Texp(ij |h(r)|sin/3(r)5i + |h(i)|cos/3(r)(r)5 2 )dr), (38) 

= Texp(i f a(T)S 1 +b(T)S 2 )dTy (39) 



We thus have 

N(t) = (a(t)S 1 +b(t)S 2 )N(t), (40) 
= ((n(t) x n(t)) • ei(t)Si + (n(t) x n(t)) • e 2 (i)5 2 )iV(i). (41) 

On the other hand, by using Property 1. in section 2, we have 

Ar\t) = A-'^nxn)^), (42) 
= (i-'WWnxhj.S^Ji- 1 !*), (43) 
= ((n(t) x n(t)) • ei(t)Si + (n(t) x n(t)) • e^)^" 1 (t). (44) 

Since JV(i) and ^4 _1 (i) satisfy the same differential equation with the same initial condition, they must be 
the same: N(t) = From this we deduce that U(t) = 1. The factorization therefore recovers the 

sudden approximation as a special case. 

The explicit form of the angular velocity of N(t) also suggests that N(t) itself can be made explicit under 
certain conditions that go beyond previously studied situations. Below we list two simple classes of explicit 
solutions of N(t). 

(i) The operator N(t) is explicit if f3(t) — (f>(t) = 0. This condition can be written more explicitly, by the 
definitions of /?(£) and (j)(t), as the following: 

kB(t) = • (n(t) x n(t)). (45) 

Under this condition, it is clear that N(t) becomes an ordinary exponential so that it is explicitly determined: 

N(t) = exp(i J \h(T)\S 2 dT) =exp(iS 2 l{t)), (46) 

where l(t) is the length of the path traversed by n(t) on the unit sphere. The condition can be seen as putting 
no constraint on how n(t) varies while giving B(t) explicitly for a given function n(i). Therefore, there are in- 
finitely many explicit examples that belong to this class. For instance, if n(t) = (sin 9 cos u>t, sin 9 sin cut, cos 9), 
then this condition automatically gives kB = —locos 9. 

This class is clearly outside the scope of perturbative analysis. We recall that near the adiabatic limit, the 
dynamical angle changes much more rapidly than (i(t) does so that the oscillating sine and cosine functions 
can make sure, under quite general conditions, that the small effect of |n(i)| ~ e does not accumulate 
during an adiabatic process with t G [0, 1/e]. The expression for N(t) in this class of examples goes beyond 
adiabatic perturbation theory in that even if |h(f)| is small, it may accumulate over time (because there are 
no oscillating factors in the intergrand) so that transitions among instantaneous eigenstates can be large. In 
fact, we see that in this case the transition probability among instantaneous eigenstates is characterized by 
the length of the path traversed by n(t) at time t. 

In general, n(t) can be an arbitrary smooth curve on the 2-sphere and it is also possible that B(t) changes 
sign during the time evolution. This only agrees with our assumption that the magnetic field is allowed to 
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change in a general way which includes the situation where the magnetic field B(f ) = B(t)n(t) passes through 
the degeneracy point of the Hamiltonian. In order that B(t) and n(t) change smoothly, B(t) is allowed to 
become zero and to change sign across the degeneracy point. 

(ii) The second class of explicit solutions for N(t) is obtained under the conditions that 

|n(t)|=ci, 0(t) - 4>(t) = c 2 t, (47) 

where ci and c 2 can be arbitrary constants. If these two conditions are satisfied, we have 

N = idiSt sinc 2 t + S 2 cosc 2 t)iV = i Cl e lC2tS3 S 2 e~ iC2tSs N. (48) 

Let 

N = e iC2t53 iVi. (49) 

Then we have 

JVi = (ici5 2 - *C2S 3 )JVi- (50) 
From this we see that Ni is explicitly solved and an explicit expression for N(t) is obtained: 

N(t) = e iC2tS3 e ( iciS2 ~ iC2S 3) t (51) 

This class of explicit factorizations include as a special case the well-known exact solution developed 
specifically for analyzing the phenomenon of magnetic resonance. In this case the function n(i) is simply 
n(t) = (sin#cosix>i, sin6>sincjt, cos 6) and so we have kB = — uj cos 8 — c 2 , which is an arbitrary constant. 
The traditional method for solving this special case is different from the factorization of time evolution we 
studied here. This is hardly surprising because this special solution was found long before the discovery of 
the geometric phase. Let us note that this special case also belongs to the first class if we choose c 2 = 0. 
Here, one must carefully distinguish between the condition kB = — uj cos 9 and the traditional magnetic 
resonance condition. In the traditional magnetic resonance analysis, one considers transitions among energy 
levels along the direction (0,0,1) fixed by a static magnetic field. The energy levels correspond to the 
static magnetic field only, not the total Hamiltonian, which also contains an oscillating magnetic field. In 
our notation, this static magnetic field is then —kBcosO. (Let us assume k < here.) The traditional 
magnetic resonance condition is to = —kBcosO. The difference lies in that the condition cu = —kBcosO is 
the resonance condition in the fixed direction while the condition kB = —uj cos 9 here represents an analogous 
resonance in the moving direction represented by n(t) and the energy levels along n(i) are determined by 
the whole instantaneous Hamiltonian. These two resonance conditions only agree with each other for small 
9. When 9 is close to 7r/2, the dominating effect would be transitions along the moving direction n(t). It 
may be interesting to study if the condition kB — —ui cos 9 can have experimental consequences. 

We want to point out that there are new examples that belong to this class. Now we shall give a new 
explicit example of n(t) and B(t) that satisfy Equation (47). First we note that the condition f3(t) — (j>(t) = c 2 t 
amounts to kB(t) = f3(t) — c 2 = — • ( n (^) x n(t)) — c 2 , so that we can focus on finding n(t) that satisfies 
|n(t) | = c\. Let us take n(t) = (sin At cos </?(£), sin At sin <p(£), cos A£). Then it is easily shown that 

|n(t)| 2 = A 2 + ^ 2 sin 2 At = C 2. (52) 

For A > 0, c\ > A 2 and < t < the solution for ip(t) that satisfies <p(0) =0 can be found to be 

¥>(*) = ±y c 2 -A 2 ln(tan(^ + J)). (53) 

From these solutions, kB(t) — — r^mp ■ ( n (0 x n(t)) — c 2 is also explicitly known. In this example, it is clear 
that the curve n(t) in general does not lie in a plane. 

It could be of interest to find more examples within or beyond these two simple classes listed above and 
to investigate possible experimental manifestations of the solutions obtained. 
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5 Concluding Remarks 



The main purpose of this paper is to clarify N(t) for a much studied physical model: spin in a general time 
varying magnetic field. Our motivation is to understand, in this specific example, the interplay between 
the geometric effect and the dynamical effect (as represented by the operators A(t) and D{t) respectively) 
in determining the structure of the N(t). Our main result is the explicit construction of the instantaneous 
angular velocity of N(t). This clarifies the structural information of N(t) and facilitates the study of 
transition probabilities among instantaneous eigensates for a general time evolution. It puts N(t) on an 
equal footing with the operator A(t) and D(t) in the sense that the instantaneous angular velocities of all 
these operators are now explicitly known. With this structural information of N(t), we show that N(t) 
itself can be made explicit under two classes of simple conditions, both of them go beyond the scope of 
perturbation theory. This gives new physical examples of explicit functions of the magnetic field B(£) that 
correspond to the explicit construction of the time evolution operator. 

The study of the extension of the adiabatic factorization of U(t) in the above spirit can be useful in 
other physical situations too. In a recent work, we examined the extension of the adiabatic factorization of 
the time evolution for the Landau problem with a spatially uniform electric field E(i) that has a general 
time dependence, including the situation where E(i) changes direction in time [9]. The factorization in this 
case shows that all the three factors can be determined explicitly: The non-adiabatic factor summarizes all 
non-adiabatic effects in the form of a path-ordered coherent state displacement operator while the geometric 
operator is a path-ordered magnetic translation. The path-ordered operators imply the existence of two 
numerical phase factors that are quantum mechanical in nature. And a simple expression is obtained for 
non-adiabatic transitions among different energy levels (inter Landau level transitions) for a general time 
evolution using the non-adiabatic operator. This time dependent Landau problem has been analyzed before 
from the point of view of constructing the quantum propagator [12] . Such a procedure of constructing the 
propagator however, does not seem to point to a factorization of the time evolution operator that can make 
the time evolution transparent. 

It could be of interest to study the factorization of time evolution for other physical systems, classical 
or quantum mechanical, where the same line of ideas applied here may contribute to the elucidation of the 
physical and/or mathematical aspects involved. 
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